Exciton formation in graphene bilayer 
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Exciton instability in graphene bilayer systems is studied in the case of a short-ranged Coulomb 
interaction and a finite voltage difference between the layers. Self-consistent exciton gap equations 
are derived and solved numerically and analytically under controlled approximation. We obtain 
that a critical strength of the Coulomb interaction exists for the formation of excitons. The critical 
strength depends on the amount of voltage difference between the layers and on the inter-layer 
hopping parameter. 
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I. INTRODUCTION 

Graphene, layers of two-dimensional honeycomb-array 
of carbon atoms, has attracted much interest these last 
few years due to its recent experimental accessibility 1-3 
and a wide variety of interesting properties 4-6 . Both the 
single-layer and the multi-layer graphenes are studied in- 
tensely. Much of the peculiar properties of the graphene 
layers arises from the energy spectrum near the so-called 
Dirac nodal points and the non-trivial topological struc- 
ture of the wave functions around them 7 ' 8 . 

As the engineering application of the graphene layers 
attracts increasing significance, we need to explore, ex- 
perimentally and theoretically, ways to enrich graphene's 
electrical properties and to control them. One way to 
achieve some control over the electrical properties is to 
change the number of layers and/or the bias applied 
across the layers. A recent experimental realization of 
the biased graphene bilayer is such an example 9-11 . By 
applying a gate bias across the two-layered graphene, the 
authors of Refs. 9 and 11 have observed a tunable en- 
ergy gap varying with the bias (see Fig. 1 for the bilayer 
graphene energy bands in the presence of bias). The bias 
can also potentially control the formation of excitons. 
Since the applied bias leads to the charge imbalance in 
the two layers, it is natural to suspect that the Coulomb 
attraction of the excess electrons and holes on opposite 
layers would lead to an exciton instability similar to the 
situations considered in an earlier literature 12-14 . If so, 
it will provide an additional control over the graphene as 
the formation of excitons is known to affect the electrical 
properties significantly 13,15 . 

Recent works on the exciton instability in a single- 
layer graphene arc based on the Dirac Hamiltonian 
description 16,17 . The exciton gap is derived and solved 
through a self-consistent equation similar to the one ap- 
pearing in the chiral symmetry breaking phenomenon 18 . 
It was shown that an exciton can be formed under a 
strong long-ranged particle-hole interaction 19 . Exciton 
can also be formed in a single-layer graphene through 
the mechanism of magnetic catalysis of dynamical mass 



generation as pointed out in 20. This work showed that 
the magnetic catalysis can induce exciton condensation 
even for weak particle- hole coupling 21 . These results are 
obtained in the framework of quantum electrodynamics 
QED deduced from the linear energy spectrum of the 
graphene monolayer. 

In the case of a bilayer, additional excitonic channels 
become possible as the excess electrons and holes from 
the two layers can form a "real-space" exciton. In this 
paper, we consider the possibility of an excitonic insta- 
bility in the biased graphene bilayer in the framework of 
Hartree-Fock theory. A conventional Hartree-Fock treat- 
ment had been used in the past to understand the exciton 
formation in semiconductors with success 15 . It is shown 
that the exciton can be formed if the strength of the 
Coulomb interaction U is larger than a threshold value 
U c which, for realistic graphene parameters, is compara- 
ble to the intra-layer hopping energy. The threshold U c 
is, in turn, bias dependent and can be tuned to a mini- 
mum value for an optimal bias V . Moreover, a reduction 
of the inter-layer hopping perhaps through intercalation 
is shown to greatly reduce the threshold value U c . 

In identifying excitonic channels, we consider two pos- 
sible scenarios. One is the pairing through the shortest- 
distance neighbors between the layers (a — d dimer in Fig. 
2), and the other, through the second shortest-distance 
neighbors between the layers (a — c and b — d dimers in 
Fig. 2). For each scenario we identify the threshold in- 
teraction strength U c and its dependence on the bias and 
the inter-layer hopping parameter. 

This work is divided into the following sequence. Sec- 
tion II describes the graphene bilayer and its model 
Hamiltonian including the short-range Coulomb interac- 
tion across the layers. Two excitonic channels we will 
consider in this paper are introduced. In the following 
two sections, each of these possibilities are examined in 
detail using the appropriate gap equations and their so- 
lutions. The work is summarized in section V. Some of 
the technical aspects are summarized in the Appendix. 
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FIG. 2: Graphene bilayer (Bernal stacking) 
is depicted as dashed lines. 



The a-d dimer 



FIG. 1: Energy spectrum for the graphene bilayer with t = 2.9 
eV, tx/t = 0.052 and V/t = 0.05. E^~ jt in full line and 
/t in dashed line. See text for definition of the energy 
branches labeled by Et^. 



II. FORMULATION OF THE EXCITON 
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Graphene bilayer is a two honeycomb array stacked in 
a Bernal arrangement as depicted in Fig. 2. In each layer 
the electrons can hop between nearest-neigbhour carbon 
atoms through 7r-orbitals with energy t, which is typically 
assumed at 2.9 eV. In a Bernal stacking electrons are al- 
lowed to do inter-layer hopping through the a-d dimcrs 
with the hopping energy given as 2tx with the tx/t ap- 
proximately 0.052. Here a dimer is defined as the pair of 
carbon atoms from the adjacent layers stacked along the 
c-axis. 

In writing down the Hamiltonian appropriate for the 
graphene bilayer, we denote the electron operators for 
the two sublattices in the lower layer by a% and and 
those in the upper layer by Cj and di . We assume a sym- 
metric doping due to the external bias ±V with excess 
electrons and holes on the lower and upper graphene lay- 
ers, respectively. We will be interested in the formation of 
the same-spin electron-hole exciton here, hence the spin 
degree of freedom a will be dropped. The Hamiltonian 
of the graphene bilayer in the absence of the Coulomb 
interaction reads 



H = -t \ a\bj + b)ai + c\dj + d]c 

(ij) 
i 

+Vj2(cU+4A-4a l -bib l ). (1) 



After the diagonalization (Derivation is given in Ap- 
pendix A), Eq. (1) is transformed to 



where (oik, Pk, 7fc, <5fc) now serve to define the cigenstates. 
The energy spectra depicted in Fig. 1 are the ones given 
by 



E 



±± 



2t\ 



V 2 ±2 



(3) 

The bare kinetic energy s k within the monolayer reads 
Ek = t\ J] a _i e lfc e ° |, where e^'s arc the nearest-neighbor 
vectors of the graphene monolayer: e\ = (1,0), e 2 = 
(-1/2, V3/2) and e 3 = (-1/2,-^3/2). 

The two independent nodal points i^i(2) where the 
bare electron spectrum Ek vanishes are chosen as K\ = 
' ^ 2 = m ^ ne basis (e x ,e y ) in the Bril- 

louin zone. The sum ^ a e lfc e a [ s approximately given 
by — (3/2)(k y — ik x ) near Ki and (3/2)(k y +ik x ) near K 2 - 



The bottom of the lower conduction band, E£ 



at k points where £ 
with the energy E r 



(4) 



occurs 

2 i j-2 



kirn 



txV/^W+v^ 



V 2 {V 2 + 2t 2 )/{V 2 + t\), 
The energy gap 



separating the valence and conduction bands is twice this 
value. The energy difference between the two conduc- 
tion bands or the two valence bands is \JV 2 + 2t ± — V 
at e k = 0, an d ^[4^ + t\) + W 2 i^\l(V 2 + ij) - 

Vtx/ \/V 2 + tj_ at Ek = {£k)m- These two quantities ap- 
proach t\/V and 2V, respectively, as V/tx — > oo. Gener- 
ally, the presence of both inter-layer hopping and the bias 
is essential in producing the gaps separating the various 
bands as depicted in Fig. 1. 

In describing the exciton formation, we propose to use 
the inter-layer interaction truncated to the second near- 
est neighbors as 
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Vc = U\ J~] n aii n d ,i 

i 

+U 2 ^ ( n a,in c ,i-e a + n b ,in d .i_ ect ) . (4) 

ia 

The local electronic densities are given by n a ^ = a\a,i, 
etc. The total Hamiltonian then reads H = H n +Vc- The 
U\ and U 2 terms are responsible for the exciton formation 
across the a — d dimer (nearest neighbor) and the a — c, 
b — d dimer (second nearest neighbor), respectively. 

At this point, several mean-field decoupling strategies 
present themselves. The average (a\di) might be a can- 
didate order parameter for the exciton pairing, but this 
quantity is nonzero even in the absence of any inter-layer 
interaction, provided the inter-layer tunneling t± remains 
non-zero. Only when t± = does this average become 
the exact order parameter. Nevertheless, one can use 
the "difference" (to be quantified in the next section) of 
(a\di) obtained in the presence and absence of excitons 
as the order parameter. This is the strategy we adopt to 
discuss the a — d dimer exciton formation. 

For the second-neighbor interaction, we could think of 
averages like (a\ci- ea ), and (b\di- ea ), as possible exci- 
tonic order parameters. Again, these averages are non- 
zero even in the absence of the interaction Vc- How- 
ever, since averages (ajc;_ eQ ) for a = 1,2,3 are related 
by the Z 3 symmetry, one could form linear combina- 
tions u a {a\a-e a ) which remains zero in the non- 
interacting case, but becomes a nonzero value once the 
interaction U 2 is turned on and excitons are formed. The 
appropriate linear combination is easily identified. For 
the second nearest-neighbor pairing, the excitonic order 
is directly related to the loss of Z 3 rotational symmetry 
of the lattice. 

Finally, we assume that at low energy the main mecha- 
nism of the exciton formation is due to the hybridization 
of the upper valence band (E^ ) and lower conduction 
bands (E k ), while the two outlying ones, E k + and 
E k + , remain as spectators. Accordingly the following 
reduced Hamiltonian may be used instead of Eq. (2): 



H> =Y,E k (0t0 k ~7hk), E k = E+-=-E^~. (5) 
fc 

Note that the two outlying bands are separated from the 
two inner ones by an energy difference that grows as V 
when V/t± is sufficiently large. The truncation scheme 
is expected to be valid when the bias V far exceeds the 
inter-layer tunneling energy; a situation easily realized in 
tunable gate systems 9 ' 11 . The inter-layer interaction, Eq. 
(4), will be truncated in the same subspace spanned by 
(/?fc,7fc). Such truncation greatly simplify the algebra in 
subsequent discussions. 



III. FIRST NEIGHBOR EXCITON PAIRING 

The first-neighbor interaction part reads 

V C (Ui) =U!^2 "<M»Hi = Ui^ a k+q a kA' -q d k' • ( 6 ) 
i qkk' 

According to our truncation scheme, the various oper- 
ators can be expanded in terms of (3k and 7fc operators 
corresponding to the lower conduction and upper valence 
bands, respectively. 

a k = U 12 (k)f3 k + U 13 (k)j k , 
d k = U i2 (k)/3 k + U i3 {k) lk . 

The 4x4 unitary matrix U diagonalizing the Hamiltonian 
(1) (see appendix A for a description of U) is used. The 
first-neighbor Coulomb interaction in the truncated space 
reads 

VciUr) = U 1 Y J (uUk + q)Pl +q + UUk + g) 7 t + J x 

qkk' 

{U 12 (k)p k + U 13 (k) 7k ) x 

(u: 2 (k'~q)(3l_ q + Ui 3 (k'-q) 7 l_ q ) x 

{U i2 {k')f3 k , +£/ 43 (fc')7fc')- (7) 

As our main concern is to explore the possibility of the 
excitonic order represented by nonzero (7^/?^), we will 
only keep terms from Eq. (7) involving an even number 
of (3 and 7 operators. In a Hartree-Fock approximation 
the mean-field Hamiltonian using the exciton order pa- 
rameter for the /3~7 hybridization can be written down 
as 

Vc{Ui) = -Y^Wl-lk + KllPk). (8) 

k 

The exciton gap is related to the exciton order pa- 
rameter (7^/3fe) through 

A. - f/i {^U l2 {q)UU q )(llP q )^j Ul 2 (k)U A3 {k) 
+Ui (^U: 3 (q)U i2 ( q )(^ q [3 q )^ U l3 {k)Ul 2 {k) 

-Ui {^umuim^Ip^ u* 12 (k)u 13 (k) 

-Ux Ui2{q)UUl)(liP q )^ Ul 2 (k)U i3 (k). (9) 

Combining the kinetic part and the mean-field Coulomb 
interaction Vc one obtains the full Hamiltonian 
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k k 

(10) 

The Hamiltonian (10) can be further diagonalized by the 
2x2 unitary rotation 

( fa \ _( e*«- cos0 5k e*»* shi0 5fc \( B k \ , . 

with = Afc/|Afc|, and cos2<9 5fe = E k /S k , sin 26» 5fc = 
\A k \/£ k . In terms of the eigen-operators B k , C k and 
the eigenvalue £ k = \/Ef, + |Afc| 2 , the Hamiltonian (10) 
reads H = ^2 k £ k (Bj.B k — ClC k ). The hybridization is 
given by 

By inserting expressions of the unitary matrix elements 
(Al) and the hybridization (12) in Eq. (9), one readily 
finds that the phase of the gap function is dictated in the 
manner 

A k = e~^\A k \, j* k = ^^Zy ( 13 ) 

The phase factor in the excitonic gap has a winding of 
2-7T around one Dirac point, and —2tt around the other. 
The total winding around the circumference of the full 
Brillouin zone is therefore zero. 

A general connection between the phase singularity of 
the wave function and the singularities in the order pa- 
rameters was considered in Ref. 23. This discussion can 
also be applied to excitonic order. In two dimensions, 
the topological structure discussed in Ref. 23 is defined 
as the total number of phase winding for the whole Bril- 
louin zone, which in this case is zero. In fact, the phase 
winding around K\ and Ki in Eq. (13) can be removed 
by a gauge transformation 24 



A k /t 




FIG. 3: (color online) Exciton gap obtained from equation 
(15) at T = 0. Here the parameters are : t — 2.9 eV, V/t = 
0.1, Ui/t = 8.5, t±/t = 0.052 for a 60 x 60 lattice in the 
reciprocal space spanned by k = fciRi -l-feR^- The reciprocal 
vectors Ri and R2 are denned by Ri = ^(1, y/3) and R 2 = 

£(-l,V3). 
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FIG. 4: (color online) Zero temperature exciton gap magni- 
tude |Afc| at k — (half way between K\ and K2) depending 
on the short-ranged Coulomb interaction Ui and the bias V 
applied on the layers, obtained for 30 x 30 lattice with t = 2.9 
eV, t±/t = 0.052. The green line correspond to the threshold 
Uic(V) solution of Eq. (15) for a lattice size of 500 x 500 
carbon atoms. 



P' k =h^\ i k = lk . (14) 

With this transformation, (7^ . j3' k ) becomes real, and the 
phase vanishes at K\ and K^. 

Taking out the phase, the gap equation (9) becomes 

|A fe | = ^E( 1 - sin20 2 P )(l-sin20 2fc ) 
p 

x[l + cos(20 4p -204 fe )]^tanh f^j . (15) 

The various factors are defined in Appendix A. The equa- 
tion can be solved numerically for given values of U, V 



and the inter-layer hopping parameter t± . The numerical 
solution of Eq. (15) is depicted in Fig. 3. Regarding the 
momentum dependence of the exciton gap we see that it 
increases from zero at the nodal points to reach satura- 
tion far from K\ and K2- 

Figure 4 shows the solution of Eq. (15) for a range of 
Ui and various values of the bias V. At zero temperature 
a second order phase transition of the exciton gap takes 
place with respect to the Coulomb interaction U± for each 
given bias V. This threshold value U\ c at which excitons 
begin to form is a function of V and is shown as a green 
line in Fig. 4. U\ C (V) reaches a minimal value U\ c /t ~ 
3.5 for an "optimal" choice of the bias V\ which is found 
at Vi /t ~ 1. 
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FIG. 5: Energy average {E^, } depending on the bias V . The 
average of the energy is computed by summing over the 
whole Brillouin zone with t = 2.9 eV and t±/t — 0.052. 



Interestingly, the dependence on the bias Ui c (V) ap- 
pears to be related to the behavior of the exciton gap 
Afc =0 , obtained far away from the Dirac points at k = 0, 
as shown in Fig. 4. The non-monotonic dependence 
of the gap value on V is apparent. A similar behavior 
is observed in the conduction-valence band energy gap, 
as exemplified in the Brillouin zone average (Et ) = 



J2keBZ ^k 



shown in Fig. 5. 



Using Eq. (15) we can deduce the dependence of U C {V) 
on the inter-layer parameter t±. As Fig. 6 shows, the 
threshold value decreases with t± and tends to zero as 
tj_/t — ► 0. Reducing the inter-layer hopping parame- 
ter would reduce the threshold U\ c of the short-ranged 
Coulomb interaction above which excitons can form. In- 
tercalation of layers of non-doping and insulating atoms 
between the two carbon layers would reduce significantly 
the inter-layer hopping parameter toward zero. The con- 
comitant reduction in the Coulomb interaction Ui with 
distance will be sufficiently slow compared to the expo- 
nential decay of t±, so that the regime U > U C (V) can 
be attained for a range of bias around V\ . Our analysis 
suggests that searching for ways to reduce the inter-layer 
hopping parameter experimentally would shed more light 
on the physics of exciton formation in graphene bilaycr. 

Finally, Fig. 7 depicts the behaviour of the critical 
temperature T c with respect to the Coulomb interaction 
U and for various hopping parameters t±. 



IV. SECOND NEIGHBOR EXCITON PAIRING 



Our approach in the previous section was based on the 
interaction with U 2 = 0. The average {a\di) read 




FIG. 6: Threshold U\ c depending on the bias V at zero tem- 
perature, t — 2.9 eV and various values of t±. 
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FIG. 7: Critical temperature T c depending on the ratio Ui/t 
for several values of t±/t. The critical temperature has been 
computed for 100 x 100 sublattice size and for the optimal 
value of the bias V/t = 1. T c tends to zero when t±/t — > 0. 



it*) = 5>£d fc ) 

k 

: £ (ut 2 {k)u 42 (k)(4p k ) + uuk)u^{k){ihk 

k 

J2 (Ui2(k)U4 3 (kM lk ) + UUWMUh) 



(16) 

in the scheme where the lowest and highest energy bands 
were truncated out. This average is nonzero even without 
the excitons at arbitrary temperature, and is not a good 
measure of the possible phase transition in the model. 
Instead, we relied on the fact that (P k jk) is zero unless 
the excitons exist, and used this average as a measure 
of the excitonic order and excitonic phase transition in 
the model. And indeed this order parameter vanished at 
high enough temperature and/or weak enough coupling, 
allowing us to identify the critical points, and so forth. 



6 



In this section, we search for an excitonic order param- 
eter defined in real space, which also vanishes identically 
for a non-excitonic phase. The averages (c\_ e eij) and 
(A-e bi) are given by 



(4-e a i )=YJ k -e- i ^f{k), 



k 



k 



* h f(k), (17) 



in the non-interacting limit, Eq. (2). The three unit 
vectors e a were defined earlier as the difference of the 
nearest neighbor positions in a given graphene layer. 
Here f(k) is a function whose detailed form is unim- 
portant to us. The combination e~ l ^ fc /(fc) is symmet- 
ric under the 120° rotation of the k vector, which in 
turn implies that (cjL^a.;) = (ct^Oj) = (ct^aj), and 

This observation suggests a strategy for defining an 
appropriate order parameter. First define hi a = c\_ ea a,i 

and gi a = d\_ e bi, then one can form the following linear 
combinations 





= hn 


x? ] 


= h i2 


x (3) 


= 



fH2 - - (hn +h iS ), 5- 2) = gi2 - - (ga + g&) 



(18) 



The operators x\^ and ^% have a zero average value in 
the non-excitonic phase, U2 = 0, due to the underlying Z3 
symmetry. In turn, non-zero value of one of the averages 
implies the Z% symmetry is spontaneously broken. 

The short-ranged Coulomb interaction (4) with Ui = 
and U2 7^ will render the mean-field Hamiltonian 



((a\ C i-e a )4- ea a * + h - C ) 



U 2 J2(( b l d i-*M-eJi + h- 



(19) 



The Coulomb interaction expressed in Eq. (20) is fully Z3 
symmetric (see Appendix B for the full expression of Eq. 
(20) in terms of the operator hi only). We remark that 
the second line of Eq. (20) is irrelevant for the exciton 
formation and can be dropped. 

Assuming translational invariance we can take (xl^ ) = 
(x^) and (3$ ) = (S'' 3 * 1 ), and express the interaction 
as 



V C {U 2 ) = -Y,(&kPbk+&tllPk) (21) 



where ®k expresses the exciton gap. Using the total 
Hamiltonian H = £ fe E k {fi\p k - 7*7*) - £ fc (Qk^k + 

®fc7fe^fcj' one can derive the averages (x^) and (3^) 
of the exciton order parameter from the self-consistent 
equations 



(V 



^ 2 ~(/3) 2/r.Q \ 1 i. U 

Vfc cos^(26» 2 fe)-^ tanh — - 

k 



;cos 2 (20 lfc ) Im ( 5>f 
-cos 2 (2^)Re [^/Hx^) 



(22) 



for /? = {1, 2, 3}. The averages (3^)) are related to (x {l3) ) 
by the simple relation (3^) = — (x ) f° r any (3. We 
defined 



uife.ei ^ g ifc.e 2 _|_ e ik.e 3 



3 ifc.ei _ _ _|_ g ifc.e3^ 



(23) 



and the Z3 symmetric counterparts Vj. 2 \ ¥>j. 3 \ ^j?^; and 
<^j. accordingly. The energy of the quasi-particles reads 
£ & = |0fc| 2 and the exciton gap is given by 



In terms of the new operator x\ and ^ /3 ' ) just defined, 
it can be recast in the form 



Vc(U2)^-- 9 U2Y / {T,(^)x 



t/ 2 



COs(26> 2 /c) 



cos(2(9 lfe ) x 4i Im I ^^(x ) 

/3 



■cos(2Mx4Re ( V) 

/3 



(24) 



+ (/lie fe) + ft-C 



^20) The system of self-consistent equations (22) admits an 
ensemble of solutions all obeying ^2g{x ) — 0. As it 
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turns out, the numerical solution always follows the con- 
dition that two of the |x^l' s are the same and differ- 
ent from the third. Furthermore, the phases of the two 
equal-amplitude bonds can be made equal through phase 
re-definition of the operators, and we can choose, for in- 
stance, x = X 7^ X without loss of generality. The 
other choices are related by Z% permutation. 

We will now exclusively consider the configuration 
( X (D) ± ( x (2)) = (x (3)) wherCi duc t0 j2 pX W = 0j 

the following relation holds: 

(X (1) ) = -2( X ( 2 >)=-2( X ( 3 >). (25) 

Introducing relation (25) into Eq. (22), one gets a single 
self-consistent equation of the exciton instability 




(x<"> = D i 



k 



1 



cos"(20 2fe )— tanh^ 



icos 2 (26» u .) 



Im 



\29 ik ) Re 



(2) 



(2) 



(3) 



(3) 



(x (1) 



(26) 



Solution of this can be used to generate the exciton gap 
6 fe using Eq. (24). 

Figure 8 represents the amplitude of the exciton gap 
|0fc| over the whole Brillouin zone of the graphene bi- 
layer. The exciton gap vanishes at the Dirac nodal points 
K\ and Ki as well as for a wave vector k = |(Ri + R2) 
where Ri = ^(1, s/Z) and R 2 = ^(-1, The van- 
ishing of the exciton amplitude at the point ^(Ri + R2) 
marks the breaking of the Z 3 symmetry. 

With Eq. (26), one can derive the threshold Coulomb 
interaction strength which reads 



FIG. 8: The exciton gap amplitude |6fe|/t is plotted over the 
whole Brillouin zone of the graphene bilayer. The parameters 
are t±/t = 0.052, V/t = 1. The Coulomb interaction U% = 
and we have chosen an arbitrary value XJ-z/t = 3. For a 50 x 50 
lattice in the reciprocal space spanned by k — feiRi + A^R^- 
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FIG. 9: (color online) Contour of the amplitude of the exciton 
average (x^ 1 " 1 } m the configuration for which (x' 1 ') 7^ (x^) 7^ 
(X (3> ). Here t = 2.9 eV, t±/t = 0.052, V/t = 1 and we used a 
sublattice of 30 x 30 carbon atoms. 



(27) 



Figure 10 shows the variation of U C 2 with V/t for var- 
ious values of t±/t. The similarity of this plot to Fig. 6 
is obvious. As for the case treating the Coulomb inter- 
action on dimcr a — d we see that there is an optimal 
value Vzo/t — 1 for which the threshold Uic is minimal. 
Moreover as the inter-layer parameter t± is decreased 
(by intercalation of insulating and non-doping atoms) the 
Coulomb threshold decreases. 



The critical temperature T c follows from 



l = -^^cos 2 (20 2fc )-Ltanh 



E k 



E k 2k B T c 
^) Re -I 



- cos*(20 lfc )Im ( ) Im ( ^ - - ( ^ + ^ 



„(i) 



1 



.(2) , ,„( 3 ) 



(28) 
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FIG. 10: Coulomb threshold U^c depending on the bias V for 
various inter-layer hopping paramters t± . We used t = 2.9 eV 
and sublattices of 50 x 50 to 1200 x 1200 carbon atoms. 




FIG. 11: Critical temperature for various Coulomb interac- 
tion and different inter-layer hopping parameter t± as given 
in Eq. Here t = 2.9 eV, t±/t = 0.052, V/t = 1 and we used a 
sublattice with 50 x 50 carbon atoms. Note the similarity to 
T c plot in Fig. 7. 



Figure 11 depicts the variation of the critical temperature 
with respect to the Coulomb interaction and for various 
values of the inter-layer hopping parameter t±. As for 
the dimer a — d Coulomb interaction, exciton are formed 
at higher temperatures for smaller tj_/t. 

The behaviour observed in this section are in good 
agreement with the behaviour of the critical tempera- 
ture T c and the Coulomb threshold U\ c observed in the 
case treating the Coulomb interaction U\ on dimer a — d. 
However the Coulomb threshold U-2 C is smaller than the 
threshold U cl . For t = 2.9 cV, t ± /t = 0.052, at the op- 
timal value of the bias V /t ~ 1, one gets U\ c /t ~ 3.5, 
compared to Ui c jt ~ 1.5. 



V. CONCLUSIONS 

The graphene bilayer system was considered with a 
short-ranged Coulomb interaction acting between the 
nearest and next-nearest carbon sites in a Bernal stacking 
scheme of two carbon layers. The short-ranged Coulomb 
interaction was introduced for both nearest (U\ for a — d 
dimer) and second-nearest (U2 for a — c and b — d dimers) 
neighbors between the two layers. 

For a given bias V , or electron- hole imbalance between 
the layers, a critical Coulomb interaction strength exists 
above which the excitons form. For the first-neighbor, 
a — d dimer interaction, the critical strength is U c /t ~ 3.5 
for a bias V/t ~ 1 . The threshold becomes smaller in the 
case of only the second neighbour Coulomb interaction 
and approximately equal to U2c.lt — 1.5 at V/t ~ 1. 
Hence, doping by equal and opposite charges of the bi- 
layer system with the voltage difference applied perpen- 
dicular to the bilayer can control the excitonic properties 
of the graphene bilayer in a non-trivial way. The optimal 
value of the bias V (which gives rise to the least threshold 
value U c ) was found to be V /t ~ 1. This non-monotonic 
dependence on the bias reflects the dependence of the 
energy gap between the conduction and valence bands 
graphene bilayer on the same quantity. 

Moreover, we showed that reducing the inter-layer hop- 
ping parameter, t±_ — ► reduces the threshold near the 
optimal bias U C (V ) to zero. We suggest that intercala- 
tion of non-doping and insulating atomic layers between 
the carbon layers could reduce significantly t± in such 
a way that the screened Coulomb interaction U obeys 
the condition U > U c (for bias around the optimal value 
V ) and excitons could form. It thus seems possible that 
the formation of the exciton gap can be controlled exper- 
imentally by both applying an electric field perpendic- 
ular to the graphene bilayer and tuning the inter-layer 
hopping 11,26 . 

The next step in the study of the exciton formation 
would lie in considering the long-ranged Coulomb inter- 
action between the two carbon layers. We conjecture 
that treating the long-range Coulomb interaction might 
reduce the threshold U c toward a reasonable value acces- 
sible by real graphene bilayer systems 9 ' 25 . 
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APPENDIX A: DIAGONALIZATION OF THE 
HAMILTONIAN 

In the momentum space the bilayer Hamiltonian (1) 
reads 
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k 

where tp^ = (a k b k c k d k ) and H k = 
( 



V 




2tj_ \ 

















2t±_ 








The unitary matrix diagonalizing the Hamiltonian is 
given by a string of matrices, 



£4 = u 0k u lk u 2k u 3k , 

/ e # fc e # fc o \ 



1 

3 



1-10 
e^ fc 
\ 1 -1 / 



=ph ^V fee <^ . 
After diagonalizing with Uo k one has UQ k H k Ua k = 



( 


V+e k 














V-e k 










tj _ e i<t>k 


-V+e k 





\ 


-tie 1 *' 







-V-e k J 



The second rotation is implemented by U\ k 

O 

sin 6*i e 



/ cos 0i —sin^ie 

cos 6»i 

^ cos6»i 

V - sin 6»i e^" 

cos 2^i = — ;== sin 2^i = 



° \ 



sin 0i e 


COS 01 / 



has 



After diagonalizing with U\ k one 

u\ k Ul k H k U 0k U lk 



I e k + X t\/\ -e-^Vtx/A \ 

t^/A -£fc + A e-^*>Vtj_/A 

e^Vtx/A e fe -A tj_ /A 

\- e »PkVt ± /\ t x /A -e fc ~A / 

The next step in the diagonaliza- 



where A = y / V 2 +t 2 J _ 
tion is affected by 



U 2k = 



COS 



cos^fc -sin6> 2 fc 

sin^fe cos^fc 

cos6> 2 fc -sin6» 2 fc 

sin^fc cos0 2 fc 



, sin 28 2k 



t 



x 



v/il+ £fe 2A 2 ' 



After diagonalizing with L^fc one has 



U lk U lk U 0k H k U 0kU lk U 2k - 

( A + C fc /A 

A-5 fe /A e-'^Vtx/A 

e l< ^Vt ± /A _A + £ fe /A 

\-e^Vt x /A 



-e-^Vtx/A \ 





where £ k = y e^A 2 The final step in the diagonal- 
ization is given by U 3k = 



( cos9 3k 

cos$4 fe — sin#4 fc e 

sin0 4fe e i</>fc cos6» 4fc 





cos 26q 



cos 29a 



A 2 +£ fe 



, sin 2B 3k : 



, sin 26a 





cos 6» 3fc 



^V 2 4+(A2^)2 



After diagonalizing with one has 

Ululul k U^H k U ak U lk U 2k U 3k 

(E++ \ 



E-- 

V Ei+J 

where = ± V £ fc + A 2 + ± 2|&. Combining the 

four unitary matrices into one, U k = Uo k Ui k U 2k U3 k , the 
cigenoperators are obtained as 



= U k 



The unitary matrix U leads to the following expres- 
sions needed in the evaluation of the exciton gap, Eq. 
(9) : 





u 12 (k)u: 3 (k) 
u i2 (k)uuk) 

U 12 (k)UUk) 



ie i0fc (1-sin 26» 2 fc)(l + cos(26»i + 20 4fe )), 

-ie i ^(l-sin26» 2fc )(l-cos(26»i + 26'4fc)), 
ie i0fc (1-sin 26 2k ) sin(26»i +26 4k ), 



- e^x ( 1 - sin 26 2k ) sin(2<9i + 20 4fc ) • 



(Al) 
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APPENDIX B: SECOND NEIGHBOUR 
COULOMB INTERACTION 



The second neighbour mean-field Coulomb interaction 
with Ux = can be rewritten in terms of the operator 
hi = c\di and gt = d\bi in the following form 
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